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Abstract 

We examine consequences of the density matrix approach to quantum theory in the 
context of a model spacetime containing closed timelike curves. We find that in general, 
an initially pure state will evolve in a nonlinear way to a mixed quantum state. CPT 
invariance and the implications of this nonlinearity for the statistical interpretation of 
quantum theory are discussed. 
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1. Introduction 

Recently there has been a lot of interest in the question of whether one can formulate a well 
defined quantum field theory on spacetimes containing closed timelike curves (CTCs) and 
if so, how the theory will behave. Studies of a class of spacetimes which have well behaved 
initial and final chronal regions and a compact region of CTCs in the middle have yielded 
some interesting results. If an S matrix is defined relating the initial and final quantum 
states in the asymptotic regions, then it was found [1] that free field scattering is unitary 
but for the case of interacting fields, the Feynman propagator fails to satisfy the identities 
(analogous to the fiat space Cutkosky rules) which would establish perturbative unitarity 
of the S matrix. This paper is concerned with the various attempts that have been made 
to try and recover a Copenhagen interpretation, consistent with observations made both 
before and after the region of CTCs, for interacting fields on these spacetimes. 

In the context of generalised quantum mechanics, Hartle [2] suggests that a consistent 
probabilistic interpretation can be recovered by normalising the amplitudes with a factor 
which depends on the initial state of the system. However, in doing so, he introduces 
a fundamental nonlinearity into quantum theory. The path integral approach will be 
reviewed briefiy in Section 2. 

Deutsch [4] looked at the problem from the point of view of quantum computation 
by considering the world lines of a finite number of particles, whose quantum states are 
described by density matrices. In his approach, a density matrix corresponding to a pure 
initial state can evolve into a density matrix corresponding to a mixed quantum state. 
This shows that the evolution cannot be given by S'S'^, where 5" is a unitary matrix. In 
Section 3, we compare the evolution of density matrices (given by Deutsch's rules) in a 
spacetime containing a compact region of CTCs with the standard density matrix analysis 
in a black hole geometry. For the CTC spacetime, one might guess that the evolution 
would be described by a linear superscattering operator, $, as in the black hole case but 
one of the results of this paper is to show that the density matrix evolves nonlinearly under 
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the rules that Deutsch proposes. This nonhnearity is a much worse feature than loss of 
quantum coherence and calls into question the statistical interpretation of quantum theory. 

In a recent paper, Politzer [5] applied Deutsch's general method to a specific example, 
which has received much attention due to its simplicity. In this spacetime, two flat, space- 
like, 3 dimensional disks of radius r are located in Minkowski space at t = and t = T, 
with the same spatial coordinates. To introduce CTCs, one simply identifies the upper 
surface of the bottom disk with the lower surface of the top disk. Also, the lower surface 
of the bottom disk is identified with the upper surface of the top disk, in order to avoid 
free surfaces. In his analysis, Politzer considers the world lines of only two interacting 
fermions, one inside and one outside the radius r. This model spacetime is reviewed in 
Section 4. 

In Section 5, we examine the evolution of density matrices in the Deutsch/Politzer 
example and find that, when the initial density matrix is pure, the evolution is nonlinear 
and hence superposition is lost. CPT invariance of the time machine and the implications 
of this nonhnearity for the evolution of initially mixed states are discussed in Section 6. 

2. The Sum over Histories Approach 

This section presents a review of Hartle's approach [2] to spacetimes containing CTCs. In 
generalised quantum mechanics, developed for closed systems by Gell-Mann and Hartle 
[see, for example, 3], there is still a time ordering of field operators so in the case of 
spacetimes containing CTCs, Hartle tries to circumvent this problem by use of the sum 
over histories formalism. 

As in Feynman's sum over histories prescription, one considers amplitudes for partic- 
ular 4 dimensional field configurations given by: 

Amplitude oc e^^'"^^")!/^, 

where S is the action functional. Sets of these configurations represent coarse grained his- 
tories and a decoherence functional D[a; a') measures the amount of quantum interference 
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between two histories a, a' e A, where A is an exhaustive set of histories for the cfosed 
system. 

If D{a; a') is neghgibly small for all pairs of different cc, a' e A, then the set is said 
to decohere and probabilities can be assigned to the individual histories, given by p{a) = 
D{a; a). However, if D{a; a') ^ for any a, a' G A, then that means there is nonnegligible 
quantum interference between a and a', so there is no probabilistic interpretation for that 
particular set of coarse grained histories. 

Familiar Hamiltonian quantum mechanics can be cast in this form. The dynamics 
is described by a characteristic decoherence functional, D{a.\ a') and all the fundamental 
principles of quantum mechanics i^e.g. superposition, herniiticity, positivity) can be related 
to certain properties of D{a.\a'). Other decoherence functionals are acceptable in this 
framework as long as these basic properties are respected. Therefore, by suitably changing 
the form of D{a\ ct'), one can generalise Hamiltonian quantum mechanics. 

This is exactly what Hartle does in a recent paper [2] to try and recover a reasonable 
probabilistic interpretation of quantum theory in the presence of nonunitary evolution. 
He considers a single scalar field in a general spacetime containing a compact region of 
CTCs. The chronal regions can still be foliated with spacelike hypersurfaces and on each 
surface, a Hilbert space of states can be defined. Transition amplitudes between field 
configurations on spacelike surfaces before and after the nonchronal region are given by 
a typical path integral over all 4-dimensional field configurations which match the initial 
and final conditions. Therefore, an S matrix defining scattering through the region of 
CTCs can be defined and calculated. However, the fact that the spacetime is not globally 
hyperbolic implies that the S matrix is nonunitary (the propagator no longer has the form 
of half advanced minus retarded Green functions and hence does not satisfy the unitarity 
identities [1]). 

Hartle's proposal for a consistent probabilistic interpretation is to introduce a de- 
coherence functional which reduces to that of ordinary Hamiltonian quantum mechanics 
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when there are no CTCs, but provides a generahsed theory otherwise. This is given by: 



In the above equation, X is an operator describing the nonunitary evolution through 
the nonchronal region and p is the density matrix for the initial state of the system. 

The Cct,i3 represent chains of unitary evolution operators and projection operators in 
the chronal regions before/after the region of CTCs. The initial state evolves unitarily 
between the spacelike surfaces and the projections are on to relevant subspaces of the 
Hilbert spaces defined on these surfaces. This functional does lead to consistent assignment 
of probabilities but unfortunately the normalisation factor in the denominator renders the 
theory nonlinear in the initial density matrix p. 

If one takes a quantum cosmological point of view, where the unique initial state 
might be given by the Hartle-Hawking no boundary state, for example, then this is not 
such a problem but even so it is still an undesirable feature of the theory. For example, 
if an advanced civilisation were to create a wormhole sometime in the future, one can see 
(in principle at least) that evolved information from nonlinear effects detected now might 
prevent the civilisation from ever building the wormhole in the first place! 

Anderson [6] also proposed a rule for quantum evolution in the presence of CTCs. The 
essence of this proposal involves a 'renormalisation' of the measure density for the Hilbert 
space of states defined after the acausal region. This simple idea provides a quantum 
mechanics which is causal (in the sense that probabilities in the present are not affected by 
achronal regions of spacetime in the future) and also linear in the initial density matrix. 

However, Fewster and Wells [7] found a problem with this procedure on a physical 
level. In normal quantum mechanics, we think of a physical observable as an operator which 
is self-adjoint with respect to the inner product of the Hilbert space. In the context of 
Anderson's proposal, one would expect that observables must be represented by operators 
which are self-adjoint with respect to both inner products. However, it is shown that as 





tr(XpXt) 
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a result of this reqirement, allowed operators must commute with the nonunitary part of 
the evolution, and their expectation values will only evolve according to the unitary part. 
This is a severe restriction:- for example, one can argue that momentum is unlikely to be 
in the class of allowed observables in this proposal. 

Fewster and Wells suggest an alternative way to restore unitarity based on the theory 
of unitary dilations. They replace the nonunitary evolution operator X by a unitary 
dilation of X, which maps between enlarged (and possibly indefinite) inner product spaces. 
In essence, the extra dimensions which are introduced provide a space of states for the 
particles in the CTC region which is inaccessible to outside observers and so overall, global 
unitarity is maintained. However, when observers reduce the state on the enlarged Hilbert 
space to the part they can measure, there will be loss of quantum coherence. 

3. Mappings between tensor product Hilbert spaces 

Consider a quantum system composed of two subsystems labelled 1 and 2. If there are no 
interactions between the subsystems, the total Hilbert space of states can be written as a 
tensor product: 

In this paper, we shall be interested in mappings between Hilbert spaces of this form. If 
an observer can only measure one part of the system, then all possibilities for the second 
subsystem must be summed over and in general, the resulting state relevant for the observer 
will be mixed. 

The natural context in which to describe mixed states is given by the density matrix 
formalism. A density matrix p describing a quantum system can be written as: 

P = P^b\A){B\ , 

where summation over A and B is implied, {\A)} is a basis for H and {{B\} is a basis for 
the complex conjugate Hilbert space H. The p^B can be thought of as the components of 
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a (];) tensor on Hilbert space or alternatively as a positive semi-definite matrix with unit 
trace. In a basis in which p'^B is diagonal, we can write: 

p^B = (no implicit summation) , 

where Pa is the probability for the system to be in the pure state \A). If any = 1 
with all the others zero, then the system is said to be in a pure state. Otherwise, p 
describes a mixed state {i.e. a sum of pure states weighted with probabilities Pa such that 

In the usual situation where a system is described by density matrices p_ and in 
asymptotic past and future regions, the superscattering operator $ : p_ — > p-|- provides a 
unique linear mapping 

P+ B — ^ BC P- D 

between density matrices. S^bc^ is Hermitian in each of its pairs of indices (AB,CD), 
probability is conserved, i.e. 

(uA D X D 

^ AC = Oc , 

and when there is no loss of quantum coherence {i.e. a pure state evolves to a pure state), 
it can be factorised into a product of a standard unitary S operator and its adjoint, 

$A D qA "c d 
BC = >J C>JB 

One can see how, in this approach, nonunitarity of the S matrix is reinterpreted by saying 
that the system loses quantum coherence in the evolution and transition probabilities are 
given by the superscattering operator, $, instead of [S'p. A more detailed analysis of $ in 
the context of the Euclidean approach to quantum gravity is given in [8]. We note here 
that one can really only use $ if the evolution of density matrices is linear. In this paper, 
nonlinear evolution operators will be denoted by f2 : p_ — > p_|_. 

Locally, density matrices evolve according to the Schrodinger equation: 

Up-U^^p+ , 
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where C/ is a standard unitary evolution operator. We can write this equation in a com- 
ponent notation which wiU be useful for discussing density matrices on tensor product 
Hilbert spaces: 

TT EF ( \AB jj\ CD ( \EF 

If p_ is a density operator on 'H~ ^ and the system can be decomposed as 'H~ = Hi 07^^, 
then we can write: 

p-= Pi ® p2 or 

{P-)cD = (pr)c(P2 )g- 

Even if an operator cannot be decomposed in this form (like, for example, the operator 
U above), it is still useful to think of the first (second) index in a pair as referring to 
subsystem 1 (2). A pair of indices can be thought of as a single index referring to the 
system as a whole. 

As an example, let us consider the case of a black hole spacetime. In this example, 
we can write both the initial and final Hilbert spaces as tensor products: 

n~ ® n2 and 

After the formation of the black hole, Tig^ is interpreted as the Hilbert space of states 
for the particles hidden behind the event horizon. Before the black hole forms, we can still 
define but the density operator on this space is assumed to be in the vacuum state, 
ie P2 = |0^)(0^|. The total initial density matrix will be given by 

P- = \Ai)\02){Bi\{02l 
and in component form, the total evolution equation is 

{P+)eF =U AO {P-)bO EF- 



The particle states behind the event horizon are unobservable, so we must sum over 
all possibilities to obtain the reduced superscattering operator relevant for subsystem 1, 
the region outside the horizon: 

(hC B _TT CD 7-rt BO 
^ EA — <J AO ED ■ 

To emphasize the point made earlier, the mapping provided by this ea^ is only 
defined for subsystem 1 and the final reduced density operator obtained, p^, will in general 
correspond to a mixed state, even though the full density operator describes a pure state 

on n+. 

We now apply this formalism to the evolution of density matrices on a spacetime 
containing a compact region of CTCs. Deutsch's proposed rules for the evolution consist 
of: 

(a) an assumption that if the initial quantum state outside the time machine, p^, is 
pure at t = 0, then the total quantum state can be written as a tensor product p_ = 
Pi<S>P2- The final quantum state at t = T is not, in general, a tensor product but we can 
still obtain the density matrices relevant for the individual subsystems by tracing over the 
degrees of freedom inside/outside the time machine, 

(b) an imposed boundary condition to obtain consistency around the CTC: 

/ \EF _ I +nF jt EF I \AB jj\ CD / -\F 
iP+)EH = (P2 )h = U aB {P-)CD ^ EH = (p2 )h ■ 

We wish to obtain an expression for the operator which maps initial to final density 
matrices, just as in the black hole case. Generically, the final external state will 
correspond to a mixed state and so there is a loss of quantum coherence as in the black hole 
case. However, the imposed consistency condition creates an essential difference. When 
this boundary condition is solved, one finds that the components (p^)^ are, in general, 
a nonlinear function of so this implies that the evolution operator obtained on the 

basis of the above rules is nonlinear. 
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The final state outside the region of CTCs is given by: 



( +\^ TT EF ( ^AB j-A CD 

{Pi)g ab {P-)cdU GF 

TT EF ( -\A( -\Btt] CD 

Hence, the reduced evolution operator relevant for the state outside the time machine is 

CE C _ TT EF(-\Btt\ CD 
GA -U AB KP2 )d^ GF ■ 

By inspection, one can see that Q, depends on (p^)^ which implies that the overall scat- 
tering will also be nonlinear. 

So far, the discussion has been fairly general. In a recent paper [5], Politzer applies 
Deutscli's density matrix proposal to a specific spacetime containing CTCs. In his example, 
pure initial states generically evolve to mixed states. The evolution of initially mixed states 
is treated by employing the standard statistical interpretation, i.e. by decomposing the 
state into its constituent weighted pure states, evolving each pure state separately and then 
recombining them in the appropriate way at the end. In Section 5, we focus on a specific 
example of the nonlinear evolution discussed above in the context of Politzer 's spacetime 
and consider the implications for the evolution of initially mixed states in Section 6. First 
of all, though, we shall briefiy review the Deutsch method applied to Politzer's model 
spacetime. 
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4. The Deutsch/Politzer example 

The example consists of a fermion with fixed spatial coordinate and world line defined for 
all time t, labelled subsystem 1. For t between and capital T, where T > 0, the fermion 
can interact with an identical particle at a different fixed spatial location, subsystem 2. 
For simplicity, we have assumed that system 2 does not exist ior t > T ov t < 0. 

In this region of interaction, the Hilbert space for the system is 4-dimensional, a 
suitable basis of states is given by the {| tt),| Ti),! it); I ii)} basis and the state of the 
whole system at time t can be described by a 4x4 density matrix p{t). The great advantage 
of considering fermions in this example is that all calculations are reduced to 4x4 matrix 
manipulations. So for example, by taking partial traces, we can obtain the 2x2 matrices 
describing the state of each individual particle. For the external particle, this is given by 
Pi{t) or tT2p{t), and similarly the CTC particle is described by P2(^)- 

Certain matching conditions need to be imposed. The boundary condition on p2 that 
brings the CTC into the system is given by 



where O"*" = + e and T~ = T — e (where e is small). This is just rule (b) of Section 3 which 
enforces consistency of information around the CTC. The other matching conditions are 
given by 



P2(T-) = P2(0+) 



pi(o-) 



Pi(0+) 



P2(0-) 



P2{T+) and 



Pi{T-) 



Evolution of the density matrices is given by the Schrodinger equation 



ip = [H,p] or alternatively p{t) = U{t)p{0)U{t)~^ . 
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The evolution matrix U{t) = e is unitary, the evolution is linear in p, and in the above 
basis, the Hamiltonian can be written in matrix form as: 



H = 



1 2uJo + X 7i 72 ' 

71 uq ui 

72 UJi LOq 

V 0, 



where loq is like a mass parameter, ui corresponds to a kinetic parameter, 71 and 72 rep- 
resent Yukawa 3-particle coupling parameters, and A corresponds to a 4-fermion coupling 
parameter. 

5. Evolution of the Density Matrices 

We now look at what happens when the external particle evolves towards t = and 
beyond in a pure quantum state, according to Deutsch's prescription. So we wish to use 
the boundary condition and the Schrodinger evolution to determine pi{T~^) for a given 
pi(0~). According to rule (a) of Section 3, if Pi(0~) is pure (in other words it can be 
written as pi(0~) = \ip){tp\ where is some quantum state), then at t = the most 
general 4x4 density matrix describing the system can be written in this tensor product 
form: 

p(0-) = pi(0-)® 



a 
b* 



b 

(1-a) 



The 2x2 matrix is just p2(0"'') and a and b are undetermined parameters. 
The final output will be given by 



pi(T+)=tr2 



U 



pi(o-) 



a b 
6* (1 - a) 



and the parameters a and b are determined by the constraint equations which are just 
another way of writing the boundary condition: 



tri 



U 



a b 
b* (1 - a) 



(1-a) 

In this case when the initial state is pure, the boundary condition provides enough con- 
straints to determine pi(T+) uniquely. 
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As Deutsch noted, the boundary condition can be written as 



$P2(0+) - P2(0+) , where 



$* = tri [U [pi(O-)®*] 



(* denotes the position of the operand). 

In other words, p2 is a fixed point of $, where $ is given by the above and is a hnear 
superscattering operator on the space of density matrices describing subsystem 2. In his 
paper, he proves that every operator of this form has a fixed point so there will always 
be a solution to this boundary condition and hence in his words - "...closed timelike lines 
place no retrospective constraint on the state of a quantum system." 

We now give a specific example to show that the map from initial to final density 
matrices (for subsystem 1) is in fact nonlinear as described in Section 3. We assume the 
tensor product form: 



where pi describes the initial (pure) quantum state of the external particle and p2 describes 
the CTC particle. The final 4x4 density matrix, and hence pi(T"'"), are determined by both 
pi and p2- However, as was discussed in Section 3, P2(0"'') depends nonlinearly on pi(0"'") 
and hence the output pi{T^) will also be a nonlinear function of the initial state. 
For example, suppose we write 



where 7=1 — a, c=l — a, a and (3 are fixed parameters describing the initial state and 
a and h are initially undetermined. We evolve p using the evolution matrix 



p(0+) = pi(0+)®p2(0+) 




/I 



u = 



cosO sinO 
—sinO cos9 



\ 



1/ 
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and calculate the following constraint equations: 



tri U 




a h 
h* c 



It turns out that the parameters a and h are given by 



a = 



(a — acosO + 2 | /3 p cosi 



:0) 



and 



(1 -cos^ + 4 I p P cos^) 



b = 



Psm9 (2a - 1) 
(1 - cos^) 



i.e. there is a nonlinear dependence of a,b on a and /3 so the final state pi(r+) will be a 
nonlinear function of the initial state and superposition is lost. 

6. Conclusions 

In this paper, we have examined a proposal for evolution of density matrices (describ- 
ing initially pure states) in a spacetime containing closed timelike curves. The proposal 
assumes a tensor product form for the initial quantum state and a boundary condition 
imposing consistency of information around the CTC. We found that the evolution of the 
reduced density matrix describing the state of particles outside (but interacting with) the 
CTC region is described by a nonlinear operator O. This result poses serious problems 
for the statistical interpretation of quantum theory. In the context of Politzer's example, 
this is basically the question of what to take as the initial 4x4 density matrix when the 
external particle is in a mixed quantum state at t=0. 

Recall that in the case when the external fermion was in a pure state, the single 
tensor product form was the most general total density matrix such that tr2p = Pi- If pi 
is mixed, then now the most general total density matrix p compatible with this condition 
will not be of this simple form and it will have a greater number of free parameters. As 
a consequence, the boundary condition will not provide enough constraint equations to 
determine the final output uniquely. 
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One could (and Deutsch does) argue that a product form is still justified because of 
the lack of correlation between the 2 spins at t=0. Before t=0, the particles cannot interact 
so why should they be correlated at t=0? However, Politzer noted that from the point of 
view of particle number 2 on the CTC, it has been interacting with particle 1 in its past, 
even if those interactions were not in the other particle's past. Correlations between the 2 
systems need to be taken into consideration. 

So how could we perform the evolution if we did not assume this simple product form? 
The standard interpretation of a mixed density matrix p is to decompose it as a statistical 
ensemble of pure states weighted with probabilities ie p = ^CiPi, where pi are pure density 
matrices and the probabilities Ci obey ^ q = 1, < q < 1. So if this interpretation were 
valid, the evolved mixed state could be unambiguously determined by summing the results 
for the evolved pure states which have been weighted by the initial probabilities This 
interpretation cannot be used here because the initial to final map is nonlinear. Politzer [5] 
proposed a dynamics based on a similar decomposition which retains correlations between 
the particles at t = 0. One of his results is that a pure state evolves to a mixed state 
even for external particles which do not interact with the CTC region. In this example 
it is difficult to see where quantum coherence is being lost as there is no interaction to 
carry the information. If there was a closed timelike curve on the other side of the galaxy, 
it shouldn't affect what happens here. The explanation for this counterintuitive result is 
precisely what we have found here - the nonlinear evolution ruins any sort of statistical 
interpretation. 

So where does the proposal go wrong? Solving the boundary condition creates the 
nonlinear dependence but all the condition really enforces is consistency of information 
around the CTC, which seems reasonable. It seems more likely that the assumption of 
a tensor product form as the initial state is where the problem lies. The tensor product 
rules out the possibility of correlations between the 2 systems, but it seems likely that 
correlations are indeed present. A more realistic proposal may be given based on, for 
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example, the direct sum [7] of Hilbert spaces, in which the systems are not assumed 
independent. 

Finally, we consider whether the Deutsch/Politzer time machine is time symmetric (or 
CPT invariant). If strong CPT was obeyed, then one would expect the final total density 
matrix p{T~) to be a tensor product when p2 is the fixed point of $. However, this is not 
the case as is readily verified in simple examples. 

Weak CPT (or detailed balance) does seem to be obeyed - that is, the probability of 
transition from A to B is the same as the probability of evolving from the CPT reverse of 
B to the CPT reverse of A. This can be seen in the context of the 2 spin example quite 
easily. Under the CPT operation, a spin up (down) particle becomes a spin up (down) 
antiparticle. The inner product used to calculate transition amplitudes for density matrices 
is essentially just the trace of the dot product of the matrices so weak CPT follows from 
the cyclic property of the trace. 
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